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Abstract 

We probe the stau-neutralino co-annihilation domain of the parameter space allowed by the 
current experimental bounds on the light Higgs mass, the b —> sg decay, and the amount 
of neutralino cold dark matter within the framework of minimal SUGRA models at a 500 
GeV e + e _ linear collider. The most favorable signals of SUSY are stau pair production and 
neutralino pair production where the small mass difference between the lighter stau and the 
lightest neutralino in the co-annihilation region is ~5-15 GeV and hence generates low-energy 
tau leptons in the final state. This small mass difference would be a striking signal of many 
SUGRA models. We find that a calorimeter covering down to 1° from the beams is crucial to 
reduce the two-photon background and the mass difference could be measured at a level of 10% 
with 500 ftU 1 of data where an invariant mass of two-tau jets and missing energy is used as a 
discriminator. 



1 Introduction 


Since an international electron-positron (e _ e + ) linear collider (ILC) can measure particle masses 
very accurately, there is a growing consensus that the next high energy machine to be built 
after the Large Hadron Collider (LHC) should be an ILC. Such a machine is technically feasible 
and the initial consensus is for the TESLA design [1. The siting is still under discussion. 

There has been in the past a huge amount of analysis on methods of detecting SUSY at an 
ILC. However, the minimal supergravity (mSUGRA) model [2, 3 4J, has several special aspects 
that make its predictions clearer and more directly accessible to experimental study. Hence it 
is worthwhile to examine this particular model. The existing experiments have already begun 
to restrict the SUSY parameter space significantly. Most significant of these are the amount of 
cold dark matter (CDM), the Higgs mass bound, the b —> sy branching ratio, and (possibly) 
the muon anomaly. 

The allowed parameter space, at present, have three distinct regions [5]: (i) the stau neu- 
tralino (t-Xi) co-annihilation region where \i is the lightest SUSY particle (LSP), (ii) the Xi 
having a larger Higgsino component (focus point) and (iii) the scalar Higgs (A 0 , H °) anni¬ 
hilation funnel (2M~o — M^o ^o). These three regions have been selected out by the CDM 
Xi 

constraint. (There stills exists a bulk region where none of these above properties is observed, 
but this region is now very small due to the existence of other experimental bounds.) The dis¬ 
tinction between the above regions can not be observed in the dark matter experiments where 
only the mass of the lightest SUSY particle would be obtained. However these regions can be 
observed at the ILC or the LHC where the particles will be produced directly and their masses 
will be measured. 

The three dark matter allowed regions need very precise measurements at the colliders to 
confirm which is correct. Since the ILC is suitable for making precision measurements, the 
cosmologically allowed parameter space is under a great deal of scrutiny. In this paper we 
choose to work with the f~Xi co-annihilation region. We note that many SUGRA models 
possess a co-annihilation region and if the a M anomaly maintains, it is the only allowed region 
for mSUGRA. Coannihilation is characterized by a mass difference (AM) between t and Xi 
of about 5-15 GeV. This narrow mass difference allows the f’s to co-annihilate in the early 
universe along with the y^’s in order to produce the current amount of dark matter density 
of the universe. The co-annihilation region has a large extension for up to 1-1.5 TeV, 

and can be explored at the LHC. The main difficulty, however, in probing this region is the 
small AM value. This AM needs to be measured very accurately in order to claim that the 
co-annihilation explains the dark matter content of the universe. However, the small AM value 
generates signals with very low energy tau (r) leptons and thus makes it difficult to discover this 
region at any collider due to the large size of the standard model (SM) and SUSY background 
(BG) events. It is this question for the ILC that we address in this paper. 

At an ILC, a major source of the SM backgrounds is the large two-photon events. The pre¬ 
vious studies use counting experiments to achieve their results. The discovery significance 
is calculated using N s i gna i/ \JN B q in Ref. [BJ, while in Ref. [7], the f\ mass is measured using the 
threshold method where they either scan over various center-of-mass (CM) energies or assume 
the mass of the LSP is known from the e and the p decays (to set the beam energy) in order 
to achieve the maximum sensitivity for a given f\ mass. However, as shown in Sec. 3, we study 
the scenarios where the e and p masses are too heavy to be produced at a 500 GeV machine. 
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We investigate the accuracy of measuring AM by analyzing the shapes of invariant mass dis¬ 
tributions of two r jets and unbalanced event transverse momentum (frr)- In our present work, 
we use a fixed collider energy (y/s = 500 GeV) for the mass measurement. 

We first discuss the available mSUGRA parameter space in Sec. 2, followed by an analysis 
of the signals and cross sections in Sec. 3. Monte Carlo (MC) studies on the event selection 
cuts to probe the SUSY events and the SM background are reported in Sec. 4 and the precision 
in the mass measurements in Sec. 5. We conclude in Sec. 6. 


2 mSUGRA Parameter Space 

The models of mSUGRA depends on only four parameters and one sign. These are mo (the 
universal soft breaking mass at the GUT scale Mg); mi/ 2 (the universal gaugino soft breaking 
mass at Mg); A 0 (the universal cubic soft breaking mass at Mg); tan (3 = (H 2 ) / (Hi) at the 
electroweak scale (where H 2 gives rise to u quark masses and H\ to d quark and lepton masses); 
and the sign of //, the Higgs mixing parameter in the superpotential (W^ = fiHiH 2 ). Note that 
the lightest neutralino Xi an d the gluino g are approximately related to mi/ 2 by M~.o = 0.4 m 1( / 2 
and Mg = 2.8 m^ 2 . 

The model parameters are already significantly constrained by different experimental re¬ 
sults. Most important for limiting the parameter space are: 


• The light Higgs mass bound of M h o >114 GeV from LEP [5|. Since theoretical calcula¬ 
tions of M h o still have a 2-3 GeV error, we will conservatively assume this to mean that 
(M h o) theory >111 GeV. 

• The b — > sy branching ratio . We assume here a relatively broad range (since there are 
theoretical errors in extracting the branching ratio from the data): 

1.8 x 10” 4 < B{B -> X s x) < 4.5 x 10" 4 (1) 

• In mSUGRA the Xi I s the candidate for CDM. Previous bounds from balloon flights 
(Boomerang, Maxima, Dasi, etc.) gave a relic density bound for CDM of 0.07 < D CDM h 2 < 
0.21 (where Doom is the density of dark matter relative to the critical density to close the 
universe, and h = H /100 km/sec Mpc where H is the Hubble constant). However, the 
new data from WMAP PI greatly tightens this (by a factor of four) and the 2a bound 
is now: 

0.095 < D CDM h 2 < 0.129 (2) 

• The bound on the lightest chargino mass of M~± > 104 GeV from LEP HB- 

Xi 

• The muon magnetic moment anomaly, using both /i + and fi~ data m Using 
the e + e _ data to calculate the SM leading order hadronic contribution, one gets a 2.7 a 
deviation of the SM from the experimental result MM- (Thee+e- data appears to be 
more reliable than the r decay data and Conserved Vector Current (CVC) analysis with 
CVC breaking [TSj.) Assuming the future data confirms the anomaly, the combined 
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effects of g^ — 2 and M~± > 104 GeV then only allow /i > 0 and leave only the t-\i 

Xi 

co-annihilation domain of the relic density. 

One can now qualitatively state the constraints on the parameter space produced by the 
above experimental bounds: (a) The relic density constraint produces a narrow rising band of 
allowed parameter space in the mQ-rrii /2 plane; (b) In this band, the M h o and b —> sj constraints 
produce a lower bound on m 1/2 for all tan/3 of mi/ 2 > 300 GeV, which implies M^o > 120 

GeV and M~± > 250 GeV. 

Xi 

In the following, we will analyze the case of p > 0. In order to carry out the calculations it is 
necessary to include a number of corrections to obtain results of sufficient accuracy, and we list 
some of these here: (i) two loop gauge and one loop Yukawa renormalization group equations 
(RGEs) are used from Mg to the electroweak scale, and QCD RGE below the electroweak scale 
for the light quarks; (ii) two loop and pole mass corrections are included in the calculation 
of M h 0 ; (iii) One loop corrections to Mb and M r are included J^j; (iv) large tan [3 SUSY 
corrections to b —> 57 are included El; (v) all fi~Xi co-annihilation channels are included 
in the relic density calculation El- We do not include Yukawa unification or proton decay 
constraints as these depend sensitively on post GUT physics, about which little is known. 

Figure HI illustrates the constraints on the mSUGRA parameter space for tan (3 = 10, 40 and 
50 with Aq = 0. The narrow blue band is the region now allowed by WMAP (see Eq. 0 . The 
dotted pink lines are for different Higgs masses, and the light blue region would be excluded if 
5a^ > 11 x 10 -10 . The three short solid lines indicate the x®-p cross section values. In the case 
of tan (3 = 40 they represent (from left) 0.03 xlO -6 pb, 0.002 xlO -6 pb, 0.001 xlO -6 pb and in 
the case of tan (3 = 50 they represent (from left) 0.05 xlO -6 pb, 0.004 xlO -6 pb, 0.002 xl0 ~ 6 
pb. I 11 the case of tan [3 = 10 they represent (from left) 5 xlO -9 pb and 1 xl0 ~ 9 pb. It is 
important to note that the narrowness of the allowed dark matter band is not a fine tuning. The 
lower limit of the band comes from the rapid annihilation of neutralinos in the early universe 
due to co-annihilation effects as the light fi mass, Mf , approaches the neutralino mass as 
one lowers m 0 . Thus the lower edge of the band corresponds to the lower bound of Eq. 0 
and the band is cut off sharply due to the Boltzman exponential behavior. The upper limit 
of the band, corresponding to the upper bound of Eq. 0 arises due to insufficient annihilation 
as m 0 is raised. As the WMAP data becomes more accurate, the allowed band will narrow 
even more. (Note that the slope and position of the band changes, however as A 0 is changed.) 
Thus the astronomical determination of the amount of dark matter effectively determines one 
combination of the four parameters of mSUGRA. Since the f-y 9 co-annihilation region seems 
to be experimentally most favored (including the g — 2 effect), we probe this region. Let us 
now study the available sparticles when we try to probe this co-annihilation band in a linear 
collider. 


3 Production and Signals of SUSY Particles at an ILC 

Figure [0 shows the production cross section of 0.1 fb for (black dashed), (blue 

solid), XiX® (blue dashed-dot) and chargino pair (vertical black dot) productions. We see that 
for large tan/3 the chargino pair production is almost not observable and the selectron pair 
production is unobservable. The stau pair has the largest reach in and the neutralino 
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pair has the largest reach in rriQ. We therefore focus on the stau pair and the neutralino pair 
production cross sections. The kinematical reach of the production cross sections of f+ff and 
X 1 X 2 productions for both y/s = 500 and 800 GeV are shown in the figures. We see that the 
800 GeV ILC will have a much bigger reach. We will, however, use the 500 GeV collider to 
study the signal since it seems to be the intial CM energy for ILC. 

The possible signals for fj + f f and X 1 X 2 i n niSUGRA are the following: 

e+e" -»• n+rf ^ (r + y?) + (r“x?) (3) 

e + e~ -*• x°i X 2 ^Xi + ( tti ) -> x°i + ( t + t “ x ?) ( 4 ) 

We look at the hadronic final state of taus (rUs) in order to have larger event rates. The final 
signal thus has two Th’s plus ^t- The analysis now is quite complicated since the r’s have 
low energy due to a small AM value. We need to develop appropriate event selection cuts to 
extract the signal from the SM background which is dominated by the 7 * 7 *e + e _ . 

In general, the co-annihilation region occurs for AM ~ 5-15 GeV. We choose three points for 
rTT.i /2 = 360 GeV, m 0 = 205, 210 and 220 GeV, with A 0 = 0 and tan 3 — 40 and develop our event 
selection cuts. The masses of SUSY particles in these three representative scenarios are given 
in Tabled The values of AM for these three points are 5, 10 and 19 GeV. The first selection 
we use is the electron beam polarization. Since both signals and background processes cross 
sections are affected by it, we choose appropriate polarization to increase the significance of the 
signal. The production cross sections for y/s = 500 GeV for different polarizations are given in 
Tabled The right-handed (RH) polarization V(e~) = —0.9 enhances the 7 j + f f signal, and the 
left-handed (LH) polarization, V(e~) = +0.9 enhances the X 1 X 2 signal. The SM background, 
mentioned in table, consists of uut + t~ states arising from WW, ZZ and Zvv production and 
this background becomes smaller for a right-handed electron beam. In addition to this, we also 
have two photon processes which will be described later: e + e~ —» 7 * 7 * + e + e _ —> r + r~ (or qq) 
+ e + e _ where the final state e + e _ pair are at a small angle to the beam pipe and the qq jets 
fake a t + t~ pair. This background, does not change with beam polarization and needs to be 
suppressed by appropriate cuts. 

Table 1: Masses (in GeV) of SUSY particles in three representative scenarios of AM = Mf — 
M~ 0 for mi /2 = 360 GeV, tan/3 = 40, /x > 0, and A 0 = 0. These points satisfy all the existing 

A1 

experimental bounds on rnSUGRA. The numbers were obtained using ISA JET El- 


MC Point 

(m 0 in GeV) 

M~o 

X 2 

M fi 

M~ 0 
Xi 

AM 

1 (205) 

274.2 

147.2 

142.5 

4.76 

2 (210) 

274.2 

152.0 

142.5 

9.53 

3 (220) 

274.3 

161.6 

142.6 

19.0 


The event selection cuts with the LH polarization will be optimized to enhance the X 1 X 2 
signal and the RH cuts to optimize the f^ff signal. 
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Table 2: Cross section times branching ratio (in fb), a x B(t —> Th) 2 , for SUSY and SM 4- 
fermions (4f) production in two cases of polarizations, V{e~) = —0.9 (RH) and +0.9 (LH). The 
SUSY cross sections were obtained using ISAJET na, and WPHACT [20s was used for the cross 
sections of the vvt + t~ processes. 


r(e-) 

-0.9 (RH) 

0.9 (LH) 

SM 4f 

7.84 

89.8 

SUSY point 1 

X 1 X 2 

0.41 

6.09 


Ti + ff 

28.3 

13.2 

SUSY point 2 

XlX.2 

0.40 

6.00 


r+ff 

26.6 

12.4 

SUSY point 3 

xlxl 

0.38 

5.68 


r+ff 

23.0 

10.6 


The generation of MC samples and the analysis for the signal and the background was done 
using the following programs: (1) ISAJET (TU] to generate SUSY events; (2) WPHACT J20] for SM 
backgrounds; (3) TAUOLA [23] for tau decay; (4) Events were simulated and analysed with a LC 
detector simulation 223 ■ 

3.1 Event Selection 

In order to reduce the backgrounds we require a set of event selection cuts and these cuts are 
given in Table 0 In this table: j 2 stands for second leading r jet, p vis gives the sum of visible 
momenta and 9(j 2 ,p V i S ) is the angle between them. 9j et is the angle between a r jet and the 
beam direction. The jets are reconstructed using the JADE algorithm with Y cut > 0.0025 ra 
and selected with Ej et > 3 GeV. Such a value of the Y cut parameter helps to select narrow r-like 
jets. The jet acceptance cut is required to reduce the SM background events such as 1+1+ and 
ZZ production. The acoplanarity is defined as A = 180° — A<^>(ji, j 2 ), where A0(ji,j 2 ) is 
the azimuthal angle between two r-jets. The cut on acoplanarity is very powerful in rejecting 
two photon SM backgrounds which have a huge cross section. In order to have MC samples 
of manageable size for two photon SM processes we apply a cut Amc > 30° already at the 
generator level. (In addition for these samples we apply the generator level cut on p^ rMC > 4 
GeV and require the r to be separated from the beam line by more than 35°). We also require 
no EM clusters (a) in 5.8° < 9 < 28° where the ILC detector has no tracking system and (b) in 
the angle below 5.8° with two options of a very forward calrometer (VFD). In our calculation, 
beamstrahlung and bremsstrahlung are included in the two-photon annihilation process. The 
two photon background in our analysis is similar to that discussed in Refs. EJEj. 

The number of accepted events for each class of final states for the case > 5 , 10, and 
20 GeV are summarized in Table 0] 

• The RH polarization strongly suppresses the SM background events (I+I+ etc.) and the 
neutralino events ( XiX 2 )• We also need a 1° VFD and fa > b GeV to get a clean signal for 
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Table 3: Kinematic cuts for the LH (V = 0.9) and the RH (V = —0.9) cases 


Cut Variable (s) 

LH (V{e~) = 0.9) RH (V(e~) = -0.9) 

N jet (E jet > 3 GeV) 

Th ID 

2 

1, 3 tracks; M tracks < 1.8 GeV 

Jet acceptance 

— qj et COS 9j et < 0.7 COS 6j et \ < 0.65 

-0.8 < cos 9(j 2 ,p vis ) < 0.7 j cos0 (.72, p„i S )| < 0.6 

Missing p t 

> 5 GeV 

Acoplanarity 

> 40° 

Veto on EM clusters 

or electrons 

No EM cluster in 5.8° < 9 < 28° with E > 2 GeV 

No electrons within 6 > 28° with pt > 1-5 GeV 

Very forward calorimeter (1°(2°) - 5.8°) 

No EM cluster with E > 100 GeV 


the fa A events. With no VFD there would be approximately 4,400 SM 77 background 
events swamping the SUSY signal. 

• The LH polarization allows for the detection of the X 1 X 2 signal with fa > 20 GeV without 
a VFD (as the 77 background falls to zero then), or fa > 10 GeV with a 2° VFD. However, 
both a 1° VFD and fa > 5 GeV are necessary to detect the fa fa events and to measure 
AM in the LH case. In the case of no VFD there would be ~9,300 SM 77 background 
events with fa > 5 GeV. Note that the event selection criteria in the LH polarization 
case are different from the RH case. 

Thus we find that the VFD is essential to detect SUSY in this region of parameter space. A 
lower fa increases the number of events and the significance. A 5 GeV fa cut has been found 
to be feasible at a 500 GeV 1LC. 

It should be noted that the 1° VFD is feasible for the ILC since the TESLA design (which has 
been accepted for the ILC technology) allows a VFD coverage down to 3.2 mrad (or 0.18°) [T]. 
We also note that our study is based on head-on collisions of electron and positron. However, 
it has been shown that the VFD is still able to reduce the two-photon background events even 
in the case of a beam crossing [7] . 

The fa fa cross section has the largest reach along the co-annihilation band and one would 
use this channel to measure the mass difference. This channel needs RH polarization for 
enhancement. In Figure [21 we plot the number of events accepted in our selection criteria 
with 500 fb -1 of luminosity as a function of AM for m 0 = 203-220 GeV (mi /2 = 360 GeV) in 
the RH polarization case. We see that we have more than 100 events, which will be adequate 
for the measurement of AM as discussed in Section for AM > 4.5 GeV. Figure [21 is a plot of 
the acceptance as a function of AM for m 0 = 203-220 GeV with a 1° VFD in the case of RH 
polarization. The acceptance drops rapidly as AM goes below 5 GeV. 

The event acceptance also depends on m i/ 2 as shown in Figure El This dependence arises 
because the r’s are less energetic and its angular distribution changes as the stau becomes 
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Table 4: Number of 7777 plus ]/>t events for luminosity of 500 fb -1 for points 1, 2 and 3 
corresponding to AM = 4.76, 9.5, and 19.0 GeV, respectively. All numbers except for two- 
photon backgrounds are common for different options of VFD. 


Process 

P(e-) 
= 5 

= 0.9 (LH) 

10 20 

T(e~) = 
5 

-0.9 (RH) 
10 20 

X 2 X 1 

Pt.l 

374 

342 

260 

15 

14 

11 


Pt.2 

624 

572 

425 

26 

24 

18 


Pt.3 

743 

697 

529 

29 

28 

22 

f+f- 

Pt.l 

73 

2 

0 

122 

2 

0 


Pt.2 

524 

267 

11 

786 

437 

22 


Pt.3 

946 

781 

335 

1283 

1076 

468 

SM 4f 


1745 

1626 

1240 

129 

123 

100 

SM 77 

2-5.8° VFD 

535 

7 

0 

249 

4 

0 


1-5.8° VFD 

10 

0 

0 

4 

0 

0 


heavier. We calculate the significance (a) as Af signa i/v / ABG, where the X 1 X 2 events are also 
treated as backgrounds, for a window of M eS = M(ji,j 2 , ft) (invariant mass of two r-jets and 
missing energy). For AM = 4.76 and 19 GeV, the allowed ranges for M e s are 0-54.5 GeV and 
0-183.5 GeV, respectively. The 5a reach for the f\ mass is found to be < 215 GeV (m±/ 2 < 520 
GeV) for AM = 4.76 GeV with a 1° VFD and ]/>t > 5 GeV. For AM = 19 GeV, the 5 a reach 
of the f mass at a 500 GeV ILC is < 237 GeV (mi/ 2 < 537 GeV). 

It should be noted that onr event selection cuts are optimized for a 500 GeV machine. In the 
case of an 800 GeV ILC, the cuts need to be re-optimized based on the new SUSY backgrounds 
and machine design limitations (e.g. the lower bound on needs to be raised). 


4 Measurement of Stau Neutralino Mass Difference 

The measurement of a small AM value is crucial since it would be a key evidence of the 
existence of the t~Xi co-annihilation. We propose the variable M e g as a key discriminator of 
the signal events from its background events. We first generate the high statistics MC samples 
for the SM and various SUSY events (by changing the m 0 value) and prepare the templates of 
the M eS distributions for the SM, X 1 X 21 an d U + rf events. Figure 0 (without the data points 
for 500 fb _1 ) shows examples of such templates for two mo values for a 2° VFD in the RH 
polarization case. The SM cross section is fitted by a blue line, the stau pair by a green line 
and the neutralino pair by a red line. The stau pair production peak separates from the SM 
background as AM increases. This is because for smaller AM, the two r signal appears like 
the t’s from the two photon background and consequently this region requires a VFD coverage 
down to 1°. From Figure Owe also find that the stau pair production cross section can be 
measured upto an accuracy of ± 4% for Point 2. 

Since the data of 500 fb” 1 of luminosity will be generated in the initial run for a few years, 












Table 5: Accuracy of the determination of AM for different VFDs. 




N~+~- 

Mi 

AM(“500 fb- 1 ” 

experiment) 

m 0 

AM 

(500 RV 1 ) 

2° VFD 

1° VFD 

205 

4.76 

122 

Not Determined 

4.711:8 

210 

9.53 

787 

9-5i 13 

9-511:8 

213 

12.4 

1027 

12.5113 

12.5113 

215 

14.3 

1138 

14-5113 

14-5113 


we then generate the MC samples equivalent to 500 ftT 1 of luminosity for particular AM values 
and fit them with the template functions generated for high statistics sample. The black lines 
in Figure 0 shows the fitting of the 500 fb _1 MC samples for point 2 with the templates of 
two different mo values of 210 and 211 GeV. (Other parameters are kept at the same values as 
before.) We then compare the y 2 f° r these fits. We find that the y 2 f° r these fits is minimum 
for the mo = 210 GeV case. We use the range of m 0 = 203-220 GeV and try to fit the 500 fbG 1 
MC sample for point 2 and determine the y 2 for all these different points. We plot the y 2 of 
these fits in Figure El and find that la in the y 2 corresponds to 9.5 ± 1 GeV. The true value 
of AM for the point 2 is 9.53 GeV. We repeat the same study for different fi masses i.e. for 
different AM. For lower AM(~ 5 GeV), we need to use a VFD of 1°. The accuracy of mass 
determination for two different VFDs is summarized in Table El showing the uncertainties are 
at a level of 10%, except for AM ~ 5 GeV where it is 20% and we have 100 f 1 f f events. 

Figure 13 illustrates the fi mass reach as a function of luminosity for a 5 a discovery with at 
least 100 events for AM ~ 5 GeV, where the AM would be determined to 20% or better. We 
find that 164 GeV and 205 GeV fi masses to be 5cr reach and 20% (or better) uncertainty in 
AM measurement with 500 and 2500 fb -1 at a 500 GeV 1LC. 


5 Conclusion 

We have probed the mSUGRA and SM signals in the f-y? co-annihilation region at a 500 
GeV ILC with 500 ft> of luminosity. In this region, the mass difference AM between the fi 
and the y° would typically be 5-15 GeV for a large range of mi/ 2 . This small mass difference 
would produce very low energy taus in the final state. The dominant SM background would 
be the two-photon process. With RH e~ beams our study has focused on the t + t~ production 
because it allowed us to reach large m^j values in the allowed parameter space. We proposed 
the invariant mass of two tau jets and missing energy variable, M(j!, j 2 ,^), to determine the 
mass difference and found the accuracy would be at a level of 10% using a 1° (or 2°) VFD 
except for AM = 4.76 GeV. For AM ~ 5 GeV, a 1° VFD would be crucial to suppress the 
two-photon background and the accuracy there would be about 20% with approximately 100 
signal events. We also calculated the discovery significance of this region and determined the 
5a reach in mi/ 2 for 500 fb -1 of luminosity. 
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Figure 1: Allowed region in the m 0 -mi/2 plane from the relic density constraint for tan/? = 
10, 40, and 50 with A 0 = 0 and n > 0. The narrow blue band by the WMAP data. The dotted 
pink vertical lines are different Higgs masses, and the current LEP bound produces a lower 
bound on mi/2 for low tan/3. The brick red region depicts the b —> sy constraint for tan/3 = 40 
and 50. For tan/3 = 10, the pink region depicts the Higgs mass region M^o < 114 GeV. The 
light blue region is excluded if 5a^ > 11 x 10~ 10 . (Other lines are discussed in text.) 
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Figure 2: Number of ThThX 1 X 1 events from fj + ff (solid circles) and X1X2 (solid triangles) 
production as a function of AM (for mo = 203-220 GeV at mi/2 = 360 GeV) in the RH 
polarization case. We assume 500 fb -1 of luminosity. 



Figure 3: Total event acceptance for —> t^t^x^Xi as a function of AM for mo = 203-220 

GeV (mi/2 = 360 GeV) in the RH polarization case. 
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Figure 4: Total event acceptance for TiT ± —* ThThXiX? as a function of m 1 / 2 for AM = 5.1 
GeV in the RH polarization case. 
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Figure 5: Example of fitting a MC sample containing SM and SUSY (mo = 210 GeV, m. 1/2 
= 360 GeV) events equivalent to 500 ffcT 1 to two M eS templates for SM+SUSY (m 0 = 210 or 
211 GeV, m. 1/2 = 360 GeV). A 2° VFD is assumed. The value of y 2 /n.d.f. is minimum when 
the events from the same SUSY parameter are in the 500 fk> sample. A 77 contribution (a 
narrow distribution around 20 GeV) is apparent. The fitting is similar for 1° VFD, except the 
77 contribution is substantially reduced. 
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Figure 6: X 2 -Xmin of fitting the 500 fb 1 sample for SUSY Point 2 (m 0 = 210 GeV, mi/ 2 = 
360 GeV) with the high statistics templates is plotted as a function of AM. 
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Figure 7: The T\ mass reach with AM = 5 GeV as a function of luminosity for a 5a discovery 
with at least 100 events. 
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